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DETERMINING MODES FOR THE SURFACE QUASI-GEOSTROPHIC 

EQUATION 


ALEXEY CHESKIDOV AND MIMI DAI 


Abstract. We introduce a determining wavenumber for the surface quasi-geostrophic 
(SQG) equation defined for each individual trajectory and then study its dependence on 
the force. While in the subcritical and critical cases this wavenumber has a uniform upper 
bound, it may blow up when the equation is supercritical. A bound on the determining 
wavenumber provides determining modes, which in some sense measure the number of 
degrees of freedom of the flow, or resolution needed to describe a solution to the SQG 
equation. 
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1. Introduction 

In this paper we introduce a determining wavenumber AQ{t) for the forced surface 
quasi-geostrophic (SQG) equation 

f)f) 

(I.I) dt 

u = R^9, 

on the torus = [0, L]^, where 0<a<2, z/>0, A = V—A is the Zygmund operator, 
and 

R^e = K-^{-d2e,die). 

The scalar function 0 represents the potential temperature and the vector function u repre¬ 
sents the fluid velocity. The initial data 0(0) G L^(T^) and the force / G for some 

p> assumed to have zero average. 

The wavenumber Ae(f) is defined solely based on the structure of the equation, but 
not on the force, regularity properties, or any known bounds on the solution. We prove 
that if two complete weak solutions 0i, 02 G A°°((— oo, cxd); I?) (i.e., lying on the global 
attractor) coincide on frequencies bellow -nvSiyi{Ag^^ then 0i = 02. While in the 
subcritical and critical cases this wavenumber has uniform upper bounds, it may blow 
up when the equation is supercritical. A bound on Ag immediately provides determining 
modes, which in some sense measure the number of degrees of freedom of the flow, or 
resolution needed to describe a solution to the SQG equation. 

The first result of finite dimensionality of a flow was obtained by Foias and Prodi for 
the 2D Navier-Stokes equations (NSE) in ll24l . where it was shown that low modes control 
high modes asymptotically as time goes to infinity. Then an explicit estimate on the number 
of determining modes was obtained by Foias, Manley, Temam, and Treve in ||23, and 
improved by Jones and Titi in A related result, the finite dimensionality of the global 
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attractor of the 2D NSE, was first proved by Foias and Temam in (see Constantin, 
Foias, and Temam ifT^ for the best available bound). See also ifTTl 1211 l22l |2^ l27l and 
references therein for more results in this direction. 

Equation (II.lb with a = 1 describes the evolution of the surface temperature field in a 
rapidly rotating and stably stratified fluid with potential velocity ifTSll . Being applicable in 
atmosphere and oceanography, this model is also very interesting from the mathematical 
point of view. Indeed, the behavior of solutions to (II.lb with k = 0 in 2D and the behavior 
of potentially singular solutions to the Euler’s equation in 3D have been found similar 
both analytically and numerically (see 0 [T3] [iTl [3^ and the references therein). Since 
L°“, the highest controlled norm, is critical when a = 1, equation (O) is referred as 
supercritical, critical and subcritical SQG for 0<a<l, q; = 1 and a > 1 respectively. 
The global regularity problem of the critical SQG equation has been very challenging due 
to the balance of the nonlinear term and the dissipative term in (II.lb . This problem is 
resolved now, with several different proofs and their adaptations to the case of a smooth 
force available ll2l fBl fTbl l28l [3111^ . 

The long time behavior of solutions to the critical SQG equations have been studied in 
E [To] [T4| [T5| |201 [34l I^. The first result on the existence of an attractor was obtained 
recently by Constantin, Tarfulea, and Vicol in ifTsll . where the authors studied the long 
time dynamics of regular solutions of the forced critical SQG using the nonlinear maximal 
principle lITbl . With the assumption that the force / G L°°(T^) n iT^(T^) and the initial 
data in the authors proved the existence of a compact attractor, which is a global 

attractor in the classical sense in iT'* for s € (1, 3/2), and it attracts all the points (but 
not bounded sets) in H^. Moreover, the authors proved that the attractor has a finite box¬ 
counting dimension. 

Later, Cheskidov and Dai Jb) proved that the critical SQG equation (II.lb with a = 1 
possesses a global attractor in T^(T^), provided the force / is solely in L'^ for p > 2. As 
the first step, it is established that for any initial data in a weak (viscosity) solution is 
bounded in on any interval [fg, oo), to > 0. The main tool is an application of the De 
Giorgi iteration method to the forced critical SQG as it was done by Caffarelli and Vasseur 
in 121 in the unforced case. This is the only part that requires the force to be in for some 
p > 2. Second, in the spirit of Q, the Littlewood-Paley decomposition technique is used 
to show that bounded weak solutions have zero energy flux and hence satisfy the energy 
equality. The energy equality immediately implies the continuity of weak solutions in L^. 
In the third step, an abstract framework of evolutionary systems introduced by Cheskidov 
and Foias Gl was followed to show the existence of a weak global attractor. Finally, with 
all the above ingredients at hand, an abstract result established by Cheskidov in a was 
applied to prove that the weak global attractor is in fact a strongly compact strong global 
attractor. 

In a very recent paper lITOl . Constantin, Coti Zelati, and Vicol showed that the 
attractor obtained in iia is indeed a global attractor in the classical sense, i.e., it attracts 
bounded sets in H^. The main ingredient here is an estimate of a (7“ norm of a solution in 
terms of the L°° norms of the solution and the force, which was done using the Constantin- 
Vicol nonlinear maximal principle M- Since the L°° is known to be bounded thanks to 
the De Giorgi iteration method, this automatically gives an absorbing ball in (7“, which 
in turn implies the existence of absorbing balls in and and hence asymptotic 

compactness in H^. This results in the existence of the global attractor. 

In this paper we start with introducing a time-dependent determining wavenumber Ae{t) 
defined for each individual trajectory 9{t) and then study its dependence on a and /. Given 
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a weak solution 9{t) of the SQG equation, we define 

— min{Ag . \p ll^pllr ^ '^p ^ f/, and ^ ^ IIIIQQ 

p<q 

for r G T-a- Here \q = 9q = Aq9 is the Littlewood-Paley projection of 9 (see 
Section|2]i, and 


CQ;,r 


a2(r + l)2V J 

co(a - 1)^ - 2V--7j 


0 < a < 1, 
1 < a < 2, 


T 

-^Q 


(^-l,oo), 0 < a < 1, 

l<a<2. 

\ a —1 ’ cn—1 J ’ 


for some absolute adimensional constant cq. Actually, the unit for cq is [cq] = [0]/[u], but 
the SQG equation (II.Il l is written so that 9 and u have the same unit. 

The first part in the definition of A resembles the dissipation wavenumber introduced by 
Cheskidov and Shvydkoy in H) for the 3D Navier-Stokes equation, also defined in terms 
of a critical norm, but based, i.e., the smallest one. In fS] it was shown that in some 
sense the linear term is dominant above that wavenumber. More precisely, it is enough to 
control a weak solution of the 3D Navier-Stokes equations in the inertial range, i.e., bellow 
the dissipation wavenumber, in order to ensure regularity. The dissipation wavenumber 
was also adapted to the supercritical SQG by Dai in ifT^ . where the smallest critical norm 
was used as well. 

Clearly, the determining wavenumber is much more restrictive than the dissipation 
wavenumber. First, a larger critical norm appears in the first condition of the definition 
of Ag^r- Second, Ag^r not only controls high modes, but also low modes, as can be seen 
in the second condition. From the mathematical point of view, this is due to the fact that 
there are more terms to control and less cancellations in this setting. 

In the first part of the paper we show that A is indeed a determining wavenumber. 


Theorem 1.1. Let a G (0,2) and 9i(t) and 02(f) be weak solutions of the SQG equa¬ 
tion (fO . Let A{t) = max{Ae-^^rit), Ag^^rit)} for some r G la- If 

(1-2) 01 (f) <yl(t) = 02(f) <yl(t), Vf > 0, 

then 

lim ||0i(f) - 02(f)||BO, = 0, 

t—^OO 

where I = Qf(r -f l)/2 when a G (0,1], and I = 2al{a — 1) when a G (1, 2). 

Moreover, if 9i(t) and 02(f) are two complete (ancient) bounded in Lf viscosity solu¬ 
tions, i.e., 01,02 G L°°{{—oo,oo);L'^), and 

(1-3) 01 (f) <yl(t) = 02(f) <yl(t), yt<T, 

for some T G (— 00 , 00 ], then 

6»i(Q= 02(f), Vf<r. 

Note that the second part of the theorem implies that for any solutions 0i(f), 02(f) on 
the attractor A, we have 0i = 62 provided (0i)<yi = (02)<yi, where 

A = {0(0) ; 0(f) is a complete bounded solution, i.e., 0 G L°°((—oo, 00 ); L^)}. 
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In il, Cheskidov and Dai proved that is a compact global attractor in the classical 
sense when a = 1. It uniformly attracts bounded sets in L^, it is the minimal closed 
attracting set, and it is the L^-omega limit of the absorbing ball . Clearly, this holds 
in the subcritical case a > 1 as well where we also have all the ingredients to apply the 
framework of evolutionary systems a. However, in the supercritical case a < 1, we only 
know the existence of a weak global attractor at this point. 

In the second part of the paper, using the De Giorgi iteration method, we extend the L°° 
estimate in 0 to the whole range a > 0. This argument requires the force / to be in 
for some p > 2!a and implies 


(1.4) V0eA 

V 

Note that this estimate holds for all a > 0, and it explains the choice of the space in 
Theorem ll.il Indeed, thanks to the Littlewood-Paley Theorem or simply the interpolation 


the Bii norm of a viscosity solution is bounded on the global attractor. Thus, one can take 
a limit as the initial time goes to —oo and show that the difference between two solutions 
that coincide below A is zero. On the other hand, the B^^ norm enjoys a better estimate 
than the V' norm. 

Using the L°° estimate (11.4b in the subcritical case a G (1,2) we show that for some r 




(a-ir 


for large enough t (or when 0 € A). Here we took p = 2 for simplicity. This gives the 
following bound on the number of determining modes N: 


N < 


ll/lb 

(a — 



In the critical case a = 1, the L°° estimate clearly is not enough to obtain a bound on A. 
However, combining it with the Holder estimate ||0(f)|| O'* ^ l|6'(0)l|oo + for some 
small h and large t, obtained by Constantin, Coti Zelati and Vicol in ifTOl . we show that 
there exists r such that 



for t large enough, provided 6(0) G Extending this estimate to solutions with 
initial data is an open problem. 


2. Preliminaries 

2.1. Notations. We denote by A < H an estimate of the form A < CB with some 
absolute constant C, and by ^ ~ S an estimate of the form CiB < A < C 2 B with some 
absolute constants Ci, C 2 - We write || • ||p = || • \\lp, and (•,•) stands for the L^-inner 
product. 
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2.2. Littlewood-Paley decomposition. The techniques presented in this paper rely strongly 
on the Littlewood-Paley decomposition, which we recall here briefly. For a more detailed 
description on this theory we refer readers to the books by Bahouri, Chemin and Danchin 
m , and Grafakos | 


Denote \ = \ for integers q. A nonnegative radial function x € (^“(K") is chosen 


such that 
(2.5) 

Let 

and 


x(6 = 


1, for 1^1 < I 

0, for Id >1. 


= x(C/2) - xiO 


‘PqiO = 


for 9 > 0 , 

X(d forg = -l. 


For a tempered distribution vector field u we define its Littlewood-Paley projection Ug in 
the following way: 


hq= 


2-7rfc-: 
6“ ^ 




T2 


Uq := AqU = ^ Uk^q(k)e'" T ^ LZ 
fceZ" 

where Uk is the fcth Fourier coefficient of u. Then 

OO 

U = ^Uq 
q=-l 

in the sense of distributions. To simplify the notation, we denote 


hq{y)u{x - y)dy, g >-1, 


Q 


U<Q 


= X] ^9’ HQ,R] = X 


= X 


q=-l 


P=Q+1 


|P-9|<1 


We will also use the Besov ; norm defined as 


Iklls?, 


E ii“i 


The following inequalities will be used throughout the paper: 

Lemma 2.1. (Bernstein’s inequality) Let n be the space dimension and r > s > 1. Then 
for all tempered distributions u, 

"( 7 - 7)1 


\Uq\\r < Xq 


’hq\U 


Lemma 2.2. Assume 2 < I < 00 and 0 < a < 2. Then 

I J UqA°‘Uq\Uq\'-~'^ dx > XqHUgWl- 

For a proof of Lemma 122] see EKm. 
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2.3. Bony’s paraproduct and commutator. Bony’s paraproduct formula will be used to 
decompose the nonlinear term. First, note that 

u ' %) — ^ ^ —2 * Vp “j~ ^ ^ Wp ■ ^Zr^<cp—2 E E Up ■ S/Vp!. 

p p p |p-p'|<i 

Due to (I2.5l l we have (p{^) = 0 when |^| < 3/4 or |^| > 2, and hence 

{fq9<q-2)>q+2 = 0, (/,g<q_2)<g-3 = 0, {fggq+l)>q +3 = 0, 

for tempered distributions / and g. Therefore, 

A,(u-Vu)= Aq{u<p -2 -Vvp) + Aq{Up ■ Vv<p- 2 ) 

g—l<p<g+2 O'—1^P^<3+2 

+ Aq(up ■'s/vp'). 

p>q-2 |p-p'|<l 
p'>q-2 

It is usually sufficient to use a weaker form of this formula: 

Aq{u ■ Vu) = Aq(u<p-2 ' VUp) + Aq{Up ■ Vu<p_2) 

|9-pI<2 I«-pI<2 

+ ^ ^ Aq{Up • VUp). 

p>q-2 

We will also use the notation of the commutator 
(2.6) [Ag, m<p_ 2 ■ V]vp := A,(u<p_2 ■ Vup) - M<p_2 ■ VAqVp. 

By definition of Aq we have 


[Ag, u<p-2 ■ Vjup = / hq{x - y) {u<p-2[y) - u<p-2{x)) Vvp{y) dy 

= / V/i,(x - y) (■u<p_ 2 (i/) - ■u<p_ 2 (a;)) Up(y) dy, 

where we used integration by parts and the fact that divu<p _2 = 0. Thus, by Young’s 
inequality, for any ?’ > 1, 


(2.7) 


||[Ag,'U<p_2 ■ VjUpIlT. ^ II V'U<p_2 II oo llt^p II r 
^ II V'U<p_2 II oo II t^p II r ■ 


/T2 


\z\\Whq{z)\ dz 


3. Proof of the first part of Theorem I1.1I 

Now we are ready to prove our first main result, which holds for all weak solutions of 
the SQG equation, even the ones that might not satisfy the energy inequality. 

Definition 3.1. A weak solution to (II. lb is a function 0 £ CwdO, T]; L^(T^)) with zero 
spatial average that satisfies (fTTT i in a distributional sense. That is, for any ip £ (T^ x 
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Theorem 3.2. Let a S {Q^2), andOiit), be weak solutions of the SQG equation il.Ik 
Let A(t) = inax{Ag-^^^r(t), Ag^^rit)} for some r G Xa. Let 


(3.8) 


0 <a<l, 

2a 

-l<a<2. 

a — 1 


If 


(3.9) dl{t)<A{t) — (l2{t)<A{t), VtG(Ti,T2), 

then 


(3.10) ||0i(i)-02(t)|lso, < VTi < io < i < T2, 

where c is an absolute constant. 

Proof Denote ui = R -^61 and U 2 = R^ 92 . Let w = 9i — 02, which satisfies the equation 

(3.11) wt + ui ■ \7w + uA°‘w + R^w • V02 = 0 


in the sense of distributions. By our assumption w<A{t) = 0 for t G (Ti, T 2 ). Recall that 


(3.12) 


T — 

-^rv — 


(I- 1 , 00 ), 0 < a < 1, 


( 2o! 
\^a-l 


-1 ’ 0-1 


1 < a < 2. 


Combining r G Xa with (13.81) one can verify that the conditions 

21 2 a 2 a „ 

2 f. I f. T —, ” + ”“>0 — 1, l-f — — — >0 

art r I 

are satisfied. These inequalities will be used throughout the proof. 

Projecting equation (13.1 11 1 onto the q — th shell, multiplying it by lwq\wq\^~‘^ , integrat¬ 
ing, adding up for all g > — 1, and applying Lemma yields 

- lk(^o)|lBO^ +Ciy f ||A“/'u;||^o^ dr < 

’ ’ to ’ 

— f I ^ f Aq{R-^W ■ V92)Wq\Wq\''~'^ dx dr 

q> — l 

— / I / Aq{ui ■ 'Vw)wq\wq\^~'^ dx dr 

•I to 1 Jk? 


q>-l' 


= f I dr + f J dr, 
■I to Jto 


(3.13) 
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for all Ti < to < t < T 2 . Using Bony’s paraproduct mentioned in Subsection l2.3l I is 
decomposed as 


I = ^ A,(i?-^W<p_2 • V(6»2)p)t 

9>-1|9-pI<2' 


I' ^ dx 


/R3 


-ly^ / Aq{R^Wp ■V{92)<p-2)Wq\Wq\’’ dx 


g>-l \q-p\<2 

-lyy yy f AqiR^wp-vi 02 )p)i 

=h + I 2 + h- 


'q\uJq 


\^-^dx 


These terms are estimated as follows. First, recall wcA{t) = 0. Let Q{t) be such that 
yl(t) = /L. Since r > /, we can choose m so that p + ^ + ^ = 1- Changing the 

order of summations and using Holder’s inequality, we infer 


yq\ \ ^q 


^-^dx 


\h\<iy2 yy / |Ag(i?^u;<p_2.v(02)p)z 

q>Q\q-p\<2''^^ 
p>Q+2 

= l yy yy \Aq{R^W<p- 2 -'^{ 02 )p)Wq\\Wq\’’~‘^ dx 


P>Q+2 \q-p\<2' 
q>Q 


<i yy a,ii( 02 )pIi. e 11^^111“' E PV'iu- 


p>Q+2 


I<?-Pl<2 

q>Q 


p' <p—2 


Then using the definition of Ag^r, Young’s inequality, Jensen’s inequality and the fact that 
WR'^WqWi < \\wq\\i, we obtain 


l/ll< 

^Cr.jT ^^ 

Y a; 

E ii«^9ii: 

\-^ yy xlr-\\R^wp.\\ 



p>Q+2 

|g-p|<2 

p' <p—2 




q>Q 


< 

CotjT ^ ^ 

E 

l^p 1; ^ E 

x;,--\\R^wp,\\i 



p>Q 

p'<p—2 


< 

CcK jT ^^ 

E a; ‘ 

ikpii;-' E 

xy\R^wp,\uxy-^f 



p>Q 

p' ^p— 

2 


< Ca,rt2l (y - 2 " E ^9 11^911*’ 

q>Q 



x;,\\R^wp,\\ix; 
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where we needed r <21/a in order to apply Jensen’s inequality at the last step. For I 2 we 
first change the order of summations and decompose it into two parts: 


^q \ I ‘"g 




E / q{R^Wp ■ V(6»2)<p-2)Wg 

9>Q|g-p|<2“'“" 

p>Q 

^ f \Aq{R^Wp ■y{92)<Q)Wq\\Wq\’-~‘^ dx 

p>Q|g-p|<2'’'“" 

q>Q 

+ Z E E / \Aq{R-^Wp-W{e2)(Q,p-2])Wq\\Wq\''~'^dx 


p>Q+2|9-p|<2- 
q>Q 


= hi + h2- 

Using Holder’s inequality, definition of Ag^r, and Young’s inequality for the first term, we 
obtain 


^21 ^^E E l|V(02)<Q||ooPVlMkglli-' 

P>Q |p-g|<2 
q>Q 

<Ca,rhJ2^Q\\^^'>^p\\l E 11^911;”^ 

p>Q |p-g|<2 

q>Q 

< Co,,rVl ^ \q\\Wq\\\. 

q>Q 

To estimate I 22 , we first use Holder’s inequality, change the order of summations, use 
Bernstein’s inequality 

^22<^ E E E l|V(02)p'||oo||i?^t^p|Mkg|ll-' 

p>Q+2 |p—g|<2 Q<p'<p—2 

q>Q 

E ii^(^2)pnioo E E IKIIE 

p'>Q p>p'+2 |p—g|<2 

q>Q 

EE E ii^9iiE> 

p'>Q p>p'+2 |p-g|<2 

q>Q 


and then definition of and Young’s inequality to infer 


h2 < Ca,rVl ^ \p' P“^t«pl|; E ll^^llE 

p'>Q p>p'+2 |p-g|<2 

q>Q 

< Ca,rVl \q ||Wg||;- 

q>Q 
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Since r > / > 2, we can choose m so that ^4-^ + j = l.To estimate we first integrate 
by parts, change the order of summations, and use Holder’s inequality: 


\h\<lYl / \Ag{R-^Wp{02)p)V{Wg\Wg\^ ^)\dx 

q>Q p>q-2 
P>Q-1 

Y [ \AqiR^WpiW^Wg\\Wgt^ dx 

q>Q p>q-2 
p>Q-l 

-^"^Y X] / \AgiR-^Wp{ 02 )p)Vwg\\'Wg\’'~‘^ dx 

p>QQ<g<p+2-^“^ 

+ 1'^ Y, f \Aq{R''~WQ{92)Q)yWg\\Wg\''~'^ dx 

Q<q<Q+2''^^ 

</^y]||i?^zi)p|M|(02)p||. E 

p>Q Q<<3<P+2 

Q<q<Q+2 


Then using definition of and Jensen’s inequality, we get 


\h\ < Ca,rt^P Y ^ ''ll-R'^'^plli 

'Y (i-l)ml 

P>Q Q<<?<p+2 

+ Ca,rVp\Q ^\\R-^Wq\\i Y 

>^q\\Wq\\Y 

Q<q<Q+2 

< Ca^rVl'^ Y. ^P X 

A, ' WwgWYxlzY^^^ 

P>Q Q<9<P+2 


+ Co,,rVp Y >'q\\Wq\\\ 


Q<q<Q+2 



X^^IKII^ 


q>Q 


were we used | + j > a — 1. Therefore, we have 


(3.14) 


|/| < Ca^rVl^ 




< c„ M 1 - 2T 


(l-2f-i) 


-I 


q>Q 


q>Q 


due to the choice of the parameters I and r as in (13.8b and (13.12b . 
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To estimate J, we start with Bony’s paraproduct formula 


J = ^ / /^q{{ui)<p-2 ■VWp)Wq\Wqf dx 

q>-l \q-p\<2-’'^'' 

E L ^q{{ui)p-Vw<p^2)Wq\Wq\^ dx 

q>-l \q-p\<2'^^'' 

- / / \{{ui)p ■yWp)Wq\Wq\’-~‘^ dx 


g>—1 p>q—2 

=Ji + J 2 + J 3 . 


Using the commutator notation (12.6K Ji can be decomposed as 


y^ (ui)<p_2 • v]t 

<?>-! |9-pI<2 ‘ 


/R3 


'p I 


\^-^dx 


E / (mi)< 9-2 • VwgWgItygl* 

-^V' y^ / ((ui)<p_2 - (ui)<q_2) ■ 

1 I 

g>-l 0-p <2 


k-p|<2 
=Jii + J 12 + J 13 , 


where we used the fact that X]|p- 5|<2 A^tUp = Wq. Notice that we have J 12 = 0, since 
div (ui)<g _2 = 0. Thanks to (12.7b . 

\\[Aq, (ui)<p_2 • V]w;p|li < ||V(ui)<p_2||oo|kp||z- 
Thus, the term Jn can be estimated as 


|dll| < I 


|[Aq, (mi)<p_2 • Vjwpl ItUgl' ^ dx 


q>Q \q-p\<2 
p>Q 


/B3 


<lY. H ||[A„(ui)<,.VK|M|u;,||i-^ 

q>Q \q-p\<2 
p>Q 

<iY. Y. iiv(ui)(Q,,]iuiitupiMitu,iii-i 

q>Q \q-p\<2 
p>Q 

+ ^E E llv(«i)<QlloolkplMk,lir 

q>Q \q-p\<2 
P>Q 

= Jill + dll2. 


For the first term we use Holder’s and Bernstein’s inequalities, 


Jill E 

E Apni(ui)pnuikpiMk,iii-' 

q>Q |g-p|<2 Q<p'<q 

P>Q 


^^E E 

E A;+"ii(ui)pnkiitupiMitu,iij-i 


q>Q |g-p|<2 Q<p'<q 
P>Q 
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and then the fact that ||wi||r < ||^i||r, definition of to get 


Jlll<Co,,rVlY^ \\Wp\\i\\wX ^ Y ^P' 

q>Q \q—p\<2 Q<p'<q 

p>Q 

<Ca,rI2lY\\'^X Y K' 

q>Q Q<p'^q 

<Ca,rVlYK\\'^<}\\\ ^P'-<1 

q>Q Q<p'^q 

< Ca,rVl Y 
q>Q 


The second term is estimated in a similar way; 


Ji12<IY XI X VIK^Op'lloolkpIMkglli ^ 

q>Q \q-p\<2 p'<Q 
P>Q 

^^X X X VII(^iVllooKlMKl||-i 

q>Q \q-p\<2p'<Q 
P>Q 

<Cc,MY X 

q>Q \q-p\<2 
P>Q 

q>Q 


To estimate J13, we start with splitting the summation 


|>/l3|<^V' V' / |((ui)<p_2 - (ui)<q_2) ■ VAqWpl Iwgl' ^ dx 

g>Q|9-p|<2“'“" 

p>Q 

<ZX X X / l(wi)p'l|VA,ri;p||u;,|'-^dj: 

q>Q \q-p\<2q-3<p'<q 
P>Q 

<^X X X / \iui)p'\\VAgWp\\wgX^ dx 

q>Q \q-p\<2 q-3<p'<Q 
P>Q 

^ ^ L \{ui)p'\\yAqWp\\Wq\^ ^ dx 

q>Q \q—p\<2 q—3<p'<q 

P>Q p'>Q 

= Jizi + Jl 32 - 


We use Holder’s inequality for the first term 


q>Q k—p|<2 q—^'^p'^Q 

p>Q 
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followed by definition of and Jensen’s inequality, 




Wn 


q>Q 


|9-pI<2 

p>Q 


< 


Ca,rVl Aq ‘ ||wg||; Ap ||Wp|l/Ap_^ 


q>Q 


|9-pI<2 

p>Q 


/ 


q>Q q>Q 

< Ca,riyl A“||W,|||. 
q>Q 


\ 




|9-pI<2 
\ p>Q 


/ 


Since r > / we can choose m so that ^ + 7 + ^-p = and estimate the second term as 


Jl32 < / P ||tU,|P P ApllWplU Y ll(^l)p'll^ 

q>Q k—p|<2 <?—3 <p^<<7 

p>Q p>Q 

< P lltUqlll"^ P Ap’^'^”"lkpl|/ Y K~^~~ 

q>Q k—p|<2 <?— 

p>Q p'>Q 

< P A5^'*'”™||ti;g|l, P A“r^"” 

9>Q 9~3<p^<g 

p'>Q 

g>Q 9~3 <p'<<j 

p'>Q 

< P A“||u;q||^ 

q>Q 


Again choosing m such that f + ^ + 


1*1 s'E E / 

9>Q|g-p|<2“'“ 

p>Q+2 


= 1 and using Holder’s inequality, we obtain 

|Ag((Mi)p • Vw<p-2)\ \Wq\’'~^ dx 


<l Y E / |Ag((Mi)p-Vu;<q)| |u;,|' ^ dx 

p>Q + 2|g-p|<2“'*'' 
q>Q 


<i Y E ii^'.iii-'E vii^^P'ii 

p>Q+2 l9-p|<2 p'<q 

q>Q 
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Now we use definition of and Jensen’s inequality to conclude that 


p>Q+2 


^XlVIkp'llr 

\q-p\<2 p'<q 

q>Q 


< 

Ca,rP'l 





q>Q 

p’<q 

< 

Cqi 

E ‘ 

-'EEK'IiaEV 



q>Q 

p'<q 

< 

Cq/ 


E^^IKII'’ 


q>Q 


where we used 1 + ^ — j >0. Finally, observe that J 3 enjoys the same estimate as due 
to the fact that IK'UijqUri < ll(^i)gl|ri for any ri G (l,oo]. Thus 


1^1 ^ Ca,rVl^ 






(3.15) 


q>Q 


< Ca.rl'P (1 - 2T- 


q>Q 


q ll'"glli 


due to (13.8b and (13.12b . 

Combining ( I3.13b -( l3.15b yields 

II^^WIIbo -IKMIIbo, < / 

for some absolute constants C and Ci. Recall that 


' a{r + 1) 


I = 


2 

2a 

a-V 


0 < a < 1, 
1 < a < 2, 


t'CK.r ^ 


Co 


a^{r + 1)2 V / 

co(a-1)" ( 1 - 2 "^-^) 


»(r+l) 


0 < a < 
1 < a < 


Hence, choosing co = 32 ^^ , we arrive at 

lk(f)||B 0 ^ - lk(fo)|lB<), < / I|A“/'w||bo^ dr,< f IM^BO^dT 

’ ’ ^ J to ’ O to 

for all Ti < to < t < T 2 . Combining it with Gronwall’s inequalities gives the desired 
result. □ 

Clearly, Theorem l3.2l imr)lies the first part of Theorem I 1.11 To prove the second part, we 
need to introduce viscosity solutions and show that the global attractor for such solutions 
is bounded in L°°. 
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4. L°° ESTIMATES 

The goal of this section is to obtain an explicit bound on viscosity solutions to (foi l 
when the force / is in for some p > 2la. 

A weak solution 6 {t) on [0, T] is called a viscosity solution if there exist sequences 
En —>■ 0 and 9n{t) satisfying 

f)f) 

-j—+ + e„A0„ =/, 

(4.16) at 

Un — R 9ny 


such that 9n ^ 9 in C'w([0,T]; L^). Standard arguments imply that for any initial data 
9o € there exists a viscosity solution 9{t) of (II. lb on [0, oo) with 0(0) = 0o (see ifTSl . 
for example). The solution 9{t) may enjoy some regularity depending on the force, but this 
is not needed for our argument. 

In the case of a = 1 and zero force, Caffarelli and Vasseur derived a level set energy 
inequality using a harmonic extension ||2l . Here we sketch a modification of the proof from 
m extended to all a > 0. 


Lemma 4.1. Let a > 0 and 9{t) be a viscosity solution to (II.lb on [0, T] with 0(0) € L^. 
Then for every X GM.it satisfies the level set energy inequality 

(4.17) \\\h{t 2 )\\l + v [ \\R'^h\\ldt<]-\\h{ti)\\l+ [ [ fhdxdt, 

^ Jti ^ Jti 

for all t 2 G [ti,T] and a.e. ti G [0, T]. ffere 9x = (9 — A)+ or 9x = {9 + A)_. 

Proof We only show a priori estimates. It is clear how to pass to the limit in (14.16b as 
e —>■ 0. Denote p{9) = (0 — A)+. Note that p is Lipschitz and 

ip'{9)p{9) = p{9). 

Multiplying the first equation of (II.lb by ip'{9)ip{9) and integrating over yields 

(4.18) 

+ 1 / / A°‘9(p{9) dx= fp{9) dx. 

Let f,gG C'°°(T^) such that g{x) = (/(x) — A)+. Then one can easily verify that 
{fix) - f{y)){gix) - g{y)) > {g{x) - g{y))^. 


Now, by Fubini’s Theorem, 



fix) - f{y) 

\x-y + Lj\^+°‘ 


gix) dy dx = 


fjy) - fix) 

|x-J/ + Ljj2+a 


giy)dydx. 


Note that (see ifTSll l 


A“ 


Cq:,7 


Y,pv. 




fix) - fiv) 


dy. 


V |x-y + Ljj2+a 
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Therefore 


/T 2 


/T 2 


= I P.V. 

^ ^a,r 

- 9 2 ^ 


f ifjx) - fiy)){gix) - g{y)) 
Jr 2 \x - y + Lj\'^+°‘ 

{gix)-giy)? 


dy dx 


■ dy dx 


t 2 jf 2 \x — y + 

= / |A^ g{x)\^ dx. 

Clearly this inequality should also holds for 0 and ip{9), giving 
/ A°‘9ip{9)dx> ( \A'^(p{9)^dx 


Thus, it follows from (14.18b that 


l_d 
2 dt 


'T 2 


ip"‘{9)dx + V ■ ) dx 


/T2 


1 


+iy [ |A 2 (^( 6 ()| dx < f fip{9)dx. 


Since the integral /.j ,2 V • [^(p^{9)u) dx = 0, this gives us the truncated energy inequality 
(14.17b . The case of ip{9) = (0 + A)_ is similar. 

□ 

Now we can use De Giorgi iteration to obtain explicit bounds on the L°° norm. For 
a = 1 this was done in ^ in the unforced case / = 0 , and similarly in || 6 l for / G L^. 
p > 2. Here we extend the proof in m to cover the whole range a > 0 . 

Lemma 4.2. Let a € (0, oo) and 9 be a viscosity solution of (II.lb on [0, c») with 9(0) € 
Lf and f G L^’(T^)/or some p G (2/q;, cxd]. Then, for every t > 0, 


(4.19) 




< 


l | 0 ( O )||2 


p+pQ —2 pg — 2 

110 ( 0 ) 11 ^”“-% 


p < oo, 


ii0(o)iir“ 


p = 00. 


Proof Consider the levels 

Afc = M(l - 2 -'=) 

for some M to be determined later, and denote the truncated function 

9k = (9 — Afc)+. 

Fix to > 0. Let Tk = fo(l ~ 2“^) and define the energy levels as: 

poo 

Uk = sup || 0 fe(f )||2 + 212 / ||A^ 0 fe(f )||2 dt. 

t>Tk JTk 

We take 9 = 9k and ti = s G (Tk-i,Tk), t 2 = t > Tk in the truncated energy inequality 
(14.17b . Then taking ti = s, t 2 = T > t, adding the two inequalities, taking limsup in T 
and then sup in t gives 


Uk<2\\9k(s)\\l 


iTk-i 


|/(a:)0fe(a:,T)| dxdr, 
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for a.a. s G {Tk-i,Tk). Taking the average in s on [Tk-i,Tk] yields 

2^+1 rca p poo p 


(4.20) 


Uk < 


to JTk^i 


9f.{s)dxds + 2 




\f{x)ek{x,t)\ dxdt. 


By an interpolation inequality between L°°{Lp') and ), we have that 

(4.21) 


||^'fe|lit+o(T2x[Tfc.oo)) - ’ 


where C is a constant independent of i/ and k. 

Note that 

6k-i>2~^M on {(a;,f) : 6»fc(a;,f) > 0 } . 

and hence 

2 ^ 0 ; 

Therefore, using the fact that Ok < Ok-i and (I4.211 i. we have 

2 fe+l poo) 


<- 


to 
2 fc+l 


/ / Ok{x,s)dxds 

JTk-i 

el_^{x, s)l{ek>o}dxds 


(4.22) 


<- 


to JTk-i 

22fca+l poo 


\0k-i\'^^°‘dxds 


toM°‘ 

^2koL-\-l 2+0. 

<rntL _ 

On the other hand, since / G with p > we obtain, forp' = (orp' = 1 

whenp = oo), 

(4.23) 

/‘OO n 


iTk-i 


\f[x)9k{t)\dxdt 


<1 


<1 


<11 


< 


< 


\9\l dx 


dt 

/ , \ ^tp' 


CAL 

POO / P ^ / / \ 

nk{2/p'-\-a—l) poo / P \ 

/ / \9k-i\^+P'‘^dx) dt 

JTk-i \Jt^ J 


2 fc( 2 /p'+Q;— 1 ) 

lip M2 /p'+“-1 

II jjl+(2-2p'+p'a)/2p' 


sup 


\9k-i\'^dx 


2 — 2p'+p'o 


iTk-i \JT^ 


| 6 »fc_i| 2 —da; 


dt 


^ vlVP/p'+°‘~^ 

Note that (2 — 2p' + p'a)l2p' > 0, since p > 2/a. Combining (14.201) . (14.22b and (14.23b 
yields 

1 

(4.24) 


Uk < 


vtoM^ 


c/rr^'^ + qi/ii 


2'=(^/"'+“-l).^l + (2-2p'+p'a)/(2p') 
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with a constant C independent of k. We claim that for a large enough M, the above 
nonlinear iteration inequality implies that Uk converges to 0 as fc —oo. Thus 9(to) < M 
for almost all x. The same argument applied to 6k = {9 + Xk)- also gives a lower bound. 

To prove the above claim (and automatically get an explicit expression for M in terms 
of to and Uq), first note that 0 < 0 almost everywhere if Uq = 0. Assume now Uq > 0. 
Denote S = {2 — 2p' + p'a)/{2p'). Note that 0 < ^ < a/2. Define 14 = VkUk with 
r]k = 2 '”^'? 7 o- We choose 


2 +p'(a - 1 ) 

m = max{4,---}, 


1 


m = 


(4.25) M = 


p'S 2Uo’ 

(4C')i2™(^+5)(217o)3 /2C2'"(i+^)|| 




UoiytoY 




2+p'(q-1) 


2+p'(a-l) 




U, 


2+p'{c_l) 

0 


Based on the choice of the parameters to, M, ijq, one can verify that 


Cpk 


22/ca+l 


l+a/2 




< A„i+“/2r/i+“/2 

— 2^-^ ^fc-i ’ 


fc > 1; 


ofe(2/p'+a-l) 1 


It follows from ( 14.241 ) that 

Vk = rjkUk <Cr]k 


vtoM^ 


9fc(2/p'+a —1) 


<i7ii+“/2r/i+“/2 i„i+'5r/i+'5 

— 2 'fc-1 ^k-l ^ 2 l^-l^k-l 

_^T/l+a/2 

for all k > 1. We also have Vq = pqUq < 1/2. Recalling that 0 < <5 < a/2, we arrive at 

k>l. 

It implies that 14 0 and hence 17^ —0 as fc —oo. The estimate ( 14.191 ) follows from 

K25^ . 

□ 


5. Global attractor and bounds on the determining wavenumber 
5.1. Global attractor. Thanks to the energy inequality, we have 


where Xq = 2^/L as before. Denote 

Bk2 = {9eL^ ■.\\9h<R2}, = 

i/Ao' 

Then for any solution 9(t) there exists time that depends only on 

9{t) G B[^2, Vf > 1^2. 


), such that 
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So the set 5^2 is an absorbing ball in L^. Moreover, there is a global attractor A C 5^2, 
A = {0(0) : 0{t) is a complete bounded trajectory, i.e., 0 G L°°((—oo, oo); L^)}- 


In il, in the critical case a = 1 , we proved that ^ is a compact global attractor in the 
classical sense. It uniformly attracts bounded sets in L^, it is the minimal closed attracting 
set, and it is the L^-omega limit of the absorbing ball B]^ 2 . This was done using the De 
Georgi iteration method to obtain continuity of solutions (which is automatically true 
in the subcritical case a > 1 ), and applying the framework of evolutionary systems in 
a. With all the ingredients at hand, the framework 0 gives the existence of the global 
attractor in the subcritical case as well. However, in the critical case a < 1, we only know 
the existence of a weak global attractor at this point. 

We also proved that the global attractor A is bounded in L°°. More precisely, let 




||0||oo < 

{0&Bl2 : ||0||oo<i?oo}, 


where 


(5.26) 



Lemma|42]implies that 


'||a-“/VI| 2 ' 


uX, 


«/2 


PQ -2 
p+pa—2 


1-a-- 

<Ao ’ 



A c Bl~ 


Moreover, for any solution 9{t) there exists time tL°° that depends only on || 0 (O)|| 2 , such 
that 

9{t) G i?Loo, Vf > fioo. 

So Bj^aa is an absorbing set. 


5.2. Proof of the second part of Theorem ll.il For a viscosity solution 9{t) on the global 
attractor, i.e., such that 9 G L°°((—oo, oo); L^), we have 

1_2 2 
<Roo‘Ri, 

for all t. Hence, hxing t in (13.10b and taking a limit as tg goes to —oo, we obtain the 
desired result. □ 

In the following two subsections we will derive explicit bounds on Ag for solutions 9 in 
the absorbing set Bl°° ■ 


5.3. The subcritical case a > 1. In this case Ag^r is a determining wavenumber for all 
r G ^j^)- When r = our estimates blow up. Nevertheless, we are able to pass 

to a limit as r —>■ =: rg and show that Ag^r < for some r < rg, where 

Ag(t) = min{Ag : Ap^ \\9p\\^<^v Vp > g, and A~° > ^ Ap||0p||oo < Cc,rv} 

P<9 


thanks to the following observation. 
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Lemma 5.1. Let a> 1. There exists a function r{M) € 

^e,r{M){t) < ^e{t), t G [Ti,T 2 ], 

provided 

II^IIl“(Tx[Ti.T2]) < M. 

Proof. We will choose r > 2al{a — 1) > 2, in which case ||0 p||t- < LM, where L is the 
size of the torus. Now since a—1 — 2/r>0we have that 

.a-l-S 


for 


\\0p 


Xp ^ 


< LM < Ca rVXp 


LM 


i-l-2/r 


=: N{r). 


\f A 0 > N (r) then we are done since the first condition in the definition of Ag^r is satisfied 
above N{r). Otherwise, for Ap G {Ae,N{r)) we have 


^ r ro 

Ap 


< N{r) 


2 _^ 

r ro 


as 


'0 ’ 


with ro = 


a- 1’ 


and 


Ap ^ ll^pllro ^ 
2 2 


One can verify that lini^_j^^- N(r) '■o = 1. Therefore 


A. 


t-a+i 


ll^pllr < 2^p 


30-«+7 


11^: 


Pl\ro 


— 2 Ap ^ ||0p||r(,Ap 

< |c„.rr7V(r)““^ 


2 _ 2 _ 

ro 


provided tq — r is small enough, which means Ap > Ag^r- This concludes the proof. □ 

We will now estimate Ag for 0 G To verify the first condition in the definition of 

Ag we note that 

II^pII^ < Aq^II^pIIoo < Ao^i?oo < kXp^ coKct- 1)2, 

2 

provided Ap > Aq ^ ■ 

For the second condition we estimate 

XI Ap||6»p||oo < 2A,i?oo < X^cov(a - if, 

p<<} 

provided A, > ( ( ■ 

Therefore, 


(5.27) Ag < max < Ag 


-1 


2i?„ 


2i?„ 


(a — l)^cory ’ \(a—l)^cor^ 

which is finite when a > 1. Clearly, when the force is large enough, the first term in ( I5.271 i 
dominates and 

2i?oo 


Ag Aq 


(a - Ifcov 
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When / G (15.26b gives 



I 


<Aj-“ 


Therefore 




( II/II2 \ 

\{a-iyv‘^) 




5.4. The critical case a = 1. In this section we consider the critical case a = 1 assuming 
that / £ L°° n and the initial data 0(0) £ H^. In this case it is known that there exists 
a global solution 6 £ (^([O, oo);H^) n Moreover, 

II^WIIc'- < l|0(io)||oo + Vf > fo + 

where h = min | , i| (see Theorem 4.2 and (4.19 ) in |[T0l). We 

will not keep track of the length L in this subsection in order not to overwhelm the esti¬ 
mates. Since / £ L°°, the radius of Bl<x, in (15.261) becomes i?oo ~ Mt. So when t is 
large enough we have ||0(f)||c'* ^ ^C'* ^ Now we will estimate the determining 

wavenumber 

= min I A, : A^|| 0 p||^ <Vp > g, and WI^pWoo < 


where r > 3 and cq is an absolute constant. Regarding the first condition, note that 


KWr < ||0p||oo < K’^RC'^ < ^^^A." 


(r + 1) 


2''P ’ 


provided 


' 2coi^ 


(r + 1)^ A h-2/r 

Rch 


As for the second condition. 


and h > -. 

r 


CqV 


Apll^plloc < 2XI-’^Rch < , 


P<Q 


(r + 1) 


2^'qr 


provided 

Therefore, 

Ag r ^ max ■ 


Since h < — and Rqh 


\q > 


2{r + 1 )" 

CoV 


Rch 


(r + l)2 \h-2/. 

, 11/11 = 


2(r + 1)" 

CqV 


Rch 


■, we obtain 


Ae,r ^ 


for some absolute (depending on L) constant c and some large enough r. 
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